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Abstract. Several of the inconsistencies plaguing the field of novae are resolved once we consider 
novae to be steady state super-Eddington objects. In particular, we show that the super-Eddington 
shell burning state is a natural consequence of the equations of stellar structure, and that the 
predicted mass loss in the super-Eddington state agrees with nova observations. We also find that 
the transition phase of novae can be naturally explained as "stagnating" winds. 



Introduction: Classical novae exhibit long duration super-Eddington luminosities while 
in their eruptive state. At least, this is the conclusion that should be reached when 
combining that the peak luminosity of novae (with Mi^d > O.SM©, Livio 1992) is always 
super-Eddington and that in all cases where the bolometric evolution was recovered 
(using UV observations), it was shown to decay slowly (e.g., Friedjung 1987, Schwarz 
et al. 2001, Shaviv 2001b). Thus, if classical novae are super-Eddington for durations 
much longer than any relevant dynamical time scale, two basic questions arise: 

• How can super-Eddington objects exist for durations much longer than their dynami- 
cal time scale ? 

• Even if such super-Eddington states exist, why do novae choose this state instead of 
following the known core-mass luminosity relation (CMLR)? 

A seemingly unrelated question, which we try to address as well, is 

• Why do the theoretical simulations consistently under predict the amount of ejected 
mass, as compared with actual observational determinations ofejecta mass? 

The notion that novae are super-Eddington contradicts the common wisdom usually 
invoked in which objects cannot shine beyond their classical Eddington limit, LEdch 
since no hydrostatic solution exists. In other words, if objects do pass LEdd, they are 
highly dynamic. They have no steady state, and a huge mass loss should occur since 
their atmospheres are then gravitationally unbound and should therefore be expelled. 
Thus, classical novae according to this picture, can pass Lgdd but only for a short 
duration corresponding to the time it takes them to dynamically stabilize after the onset 
of the thermonuclear runaway (TNR). This is indeed seen in detailed ID numerical 
simulations of nova TNRs, where novae can be super-Eddington but only for several 
thousand seconds (e.g., Starrfield, 1989). However, once they do stabilize, they are 
expected and indeed do reach in the simulations a state given by the CMLR which 
describes the Hydrogen shell burning state. Namely, we naively expect to find no steady 
state super-Eddington atmospheres. This, however, is not the case in nature. 



The Super-Eddington State: To existence of a super-Eddington state becomes natural, 

once we consider that: 

1. Atmospheres become unstable as they approach the Eddington limit. In addition to 
instabilities that operate under various special conditions (e.g., Photon bubbles in 
strong magnetic fields or the s-mode instability under special opacity laws), two 
instabilities were found to operate in Thomson scattering atmospheres (Shaviv 
2001a). Moreover, one of these instabilities does not depend on the boundary 
conditions and is therefore extremely general. It implies that all atmospheres will 
become unstable already before reaching the Eddington limit. 

2. The effective opacity relevant for the calculation of the radiative force on an inho- 
mogeneous atmosphere is not necessarily the microscopic opacity (Shaviv 1998). 
Instead, it is given by K^-^-^ = {FKy)y / {F)y. The situation is very similar to the 
Rosseland vs. Force opacity means used in non-gray atmospheres, where the inho- 
mogeneities are in frequency space as opposed to real space. For the special case 
of Thomson scattering, the effective opacity is always reduced. 

To summarize, we find that as atmospheres approach their classical Eddington limit, 
they will necessarily become inhomogeneous. These inhomogeneities will necessarily 
reduce their effective opacity such that the effective Eddington limit will not be sur- 
passed even though the luminosity can be super-classical-Eddington. This takes place in 
the external part of luminous objects, where the radiation diffusion time scale is shorter 
than the dynamical time scale of the atmosphere. Further inside the object, convection 
is necessarily excited such that the total energy flux may be super-Eddington, with the 
radiative part of it necessarily being sub-Eddington and with the convective flux carrying 
the difference. 

Nevertheless, one of the key features of super-Eddington atmospheres is that a wind 
will necessarily be accelerated. It is not the catastrophic wind naively expected from 
super-Eddington conditions, however, it is going to be significant. 

Super-Eddington Winds: The atmosphere can remain effectively sub-Eddington while 
being classically super-Eddington, only as long as the inhomogeneities comprising the 
atmosphere are optically thick. Clearly however, this assumption should break at some 
point where the density is low enough. From that radius upwards, the radiative force 
overcomes the gravitational pull and a wind is generated. 

The mass loss rate can then be obtained by identifying the sonic point of a steady state 
wind with the critical point, which is the radius where the radiative and gravitational 
forces balance each other. We then have M — AnR^Pcriticaiv sonic- Furthermore, this can 
generally be reduced to the form of 

^^W(T){L-L^,)^ (1) 

CVs 

where "M^ is a dimensionless wind "function". In principle, 'W can be calculated from 
first principles only after the nonlinear state of the inhomogeneities is fully understood. 
This however is still lacking as it requires elaborate 3D numerical simulations of the 
nonlinear steady state. Nevertheless, it can be done in several phenomenological models 
which only depend on geometrical parameters such as the average size of the inhomo- 



geneities in units of the scale height (P = d/lp), the average ratio between the surface 
area and volume of the blobs in units of the blob size (E), and the volume filling factor 
a of the dense blobs. For example, in the limit where the blobs are optically thick, one 
obtains that H' = 3S/32A/va|3(l - a)^ (Shaviv 2001b). Here v is the ratio between the 
effective and adiabatic speeds of sound. Thus, W depends only on geometrical factors. 
It does not depend explicitly on the Eddington parameter P. Moreover, typical values of 
W ~ 1-10 are obtained. 

The mass loss predicted by the super-Eddington theory (eq. 1) was compared with 
observations of super-Eddington objects which have good observational data. These 
were two novae which are not very fast and which have the best determined absolute 
bolometric evolution: FH-Ser and LMC 1988 #1. The theory was also applied to the 
Luminous Blue Variable star rj-Car. 

For the two novae, we find that the predicted mass loss rates agree with their obser- 
vations if W pa 10 ±5, which is clearly consistent with the theoretical estimate for W . 
The agreement is also with the temporal evolution of the velocities, if those are taken to 
be the primary absorption line component. 

Using W ^ 10±5, the mass loss equation can also be applied to Tj-Car, which is an 
entirely different object from novae (in mass, mass loss rate and duration, photospheric 
size etc), yet, the predicted integrated mass loss is in total agreement with the observed 
1-2M0 of ejecta, while the terminal velocity is consistently predicted as well. For more 
information on the models, see Shaviv (2001b) 

Steady State Shell Burning: The next step is to show why novae become super- 
Eddington to begin with. In particular, the nova shell burnings steady state should be 
given by the CMLR (Paczynski 1970), counter to observations. Given the contradiction, 
we look for new solutions to the stellar structure equations for systems with shell 
burning. This should describe the steady state of novae after they undergo a TNR. Unlike 
the standard derivation of the CMLR, we allow the atmospheres to be inhomogeneous. 
Namely, if the Eddington parameter is larger than a threshold (taken to be 0.85) the 
effective opacity is reduced. We still do not know the exact behavior of Kg// as a 
function of P. This could later be obtained by comparing the steady state obtained to 
the actual observations of novae, or through detailed 3D radiative hydro simulations. 
At this point however, we want to show that using a reasonable behavior of Ke//(P), a 
super-Eddington steady state does exist. 

We take Ke//(P > Pen; = 0.85) = KoPcnY /r and lV = \0 (somewhat different choices 
do not change the conclusions). We find that the CMLR obtains a super-Eddington 
branch. The main differences between the structure in this branch and the structure in 
the sub-Eddington branch are the following: (1) The super-Eddington branch has a SEW 
at the top of it, with the photosphere located in the wind. The sub-Eddington branch has 
no wind (though it could have one if a non-Thomson opacity is considered). Since the 
wind is "heavy" the actual luminosity at the photosphere could be sub-Eddington. (2) 
The super-Eddington branch has a convective layer that penetrates into the burning shell 
(most of the energy is actually released in the convection zone). In the sub-Eddington 
branch, the burning shell is all radiative. (3) The luminosity in the sub-Eddington branch 
is not a function of the core mass. It is a function of the core mass in the super-Eddington 
branch. 



Except for the very fast novae, the shell burning evolution can then be described as a 
steady state which slowly evolves due to mass loss (e.g., as is done by Kato & Hachisu 
1994, for sub-Eddington evolution and opacity driven winds). The preliminary results 
show that predictions for L{t), T{t), Vco and M^jecta agree with observations. 




FIGURE 1. The shell structure of a CNO shell burning white dwarf with a IQ^'^Mq envelope. Unlike 
the top panel, the bottom panel describes an atmosphere in which the radiative instability is taken into 
account, such that it can become porous. This allows the occurrence of the described steady state. The 
luminosity is super-Eddington. A convection zone arises and its inner boundary is well within the nuclear 
burning zone. A continuum driven wind is launched such that the photosphere is in a wind. 

A Transition Phase? An interesting situation can potentially arise if the mass loss 
predicted requires too large a luminosity to push the ejecta to infinity. This situation is 
in fact a realization of the prediction by Owocki and Gayley (1997), who first suggested 
that a hypothetically large mass loss rate would result with a stagnating wind. Since no 
steady state can obviously be obtained, they suggested that a variable state will exit. 

The first analysis performed is a ID numerical hydro- simulation, which will be 
elaborated elsewhere (Owocki & Shaviv 2002). A sample result is given in fig. 2. 
Although there are no indications that the structure should be purely radial, the limited 
(and preliminary) ID results appear to qualitatively give a typical transition phase 
behavior. Thus, we hypothesize that this is the origin of the transition phase (excluding 
of course the big dips which presumably arise from dust formation). 
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FIGURE 2. Preliminary results of ID hydro-simulations (Owocki & Shaviv 2002) where the luminosity 
of a system is slowly reduced, such that at some point the wind stagnates before reaching r^oa. The radial 
and temporal units are lOrsonic and rsonic/vs- The result resembles the behavior in the transition phase. 

Discussion: To summarize, several of the open questions pertaining to classical nova 
eruptions are resolved if we understand them as steady state super-Eddington objects. 
As such, the observed super-Eddington luminosities, the temperature and velocity evo- 
lution, as well as the mass loss rate are all naturally explained. Moreover, we have seen 



that by allowing atmospheres to become porous, not only is a steady state possible, it 
also becomes a viable solution of the shell burning stellar structure equations. 

One of the interesting results, is that the velocities predicted from the super-Eddington 
states are those observed as the primary absorption component. These are typically half 
as those observed for example in the "diffuse enhanced" or "Orion lines". Thus, the 
prediction is that the bulk of the mass outflow which arises in the super-Eddington 
steady state, should be with the primary velocity. Indeed, if one looks for example 
at the ejecta of FH Ser years after the eruption, the bulk velocity is similar to the 
primary component velocity. According to this picture, the higher speed components 
are therefore significantly less important in the mass loss budget. This however does not 
solve the question of their origin. 

One possibility is that high velocity components are the result of super-Eddington 
radiative driving of material that is already present outside the critical radius after 
steady state super-Eddington shell burning is established, but before the wind has had 
time to be established. As a result, this external material is accelerated with the full 
luminosity, which later falls down with radius as it is "consumed" to accelerate the 
primary component (which therefore reaches lower velocities). The pre-existing tenuous 
material could be a result of a shock wind from the early TNR/dynamic phase. It is hard 
to estimate from first principles the mass that should be contained in it, but an upper 
limit can be given. It cannot be larger than the amount of material that will exist in the 
steady state wind over a radius of order the sonic radius. 

Nevertheless, one can estimate from Shaviv (2000b) the velocity that the high com- 
ponent should have. Optically thin material sitting above the photosphere before the 
wind begins to accelerate will itself accelerate with the base flux of ToLedd- We can 
compare this to the acceleration in the primary component in which energy was used 
to up accelerate the "heavy" wind, thereby leaving only an energy of flux of Tobs^-edd 
at r — > oo. Therefore, at r — > oo the initially tenuous material will have a larger terminal 
kinetic energy such that: ^Mv^ jMv^ q -|- (Fq — Fofe)Lgjd, where Voo,o is the velocity 
of the steady state wind, that is, of the primary component. Taking the data for FH Ser 
on day 9 (which is the earliest for which velocity components are present) and nominal 
values for the FH Ser super-Eddington state from Shaviv (2000b), we obtain Voo ~ 1375 
km/s. This is similar to the 1300 — 1800 km/s seen for this component. Moreover, as a 
wind begins to buildup above the sonic point, initially lower material will be accelerated 
with a progressively lower flux. This will give rise to a Hubble flow in the tenuous fast 
component, in which material at larger radii has larger velocities. 

Another noteworthy point is that the integrated ejecta masses observed agree with the 
mass loss rate predicted by the super-Eddington theory. This implies that if one assumes 
the envelope masses at TNR to be roughly the observed ejected mass, then consistent 
nova eruptions are obtained with proper luminosity, temperature, mass loss rate and 
velocity evolution. The catch is that this envelope mass is typically an order of magnitude 
larger than the threshold masses required to trigger a TNR. The problem cannot be 
resolved by "dredging up" ten times as much mass as the envelope from the WD 
because the hydrogen fraction seen in the ejecta is not extremely small. Thus, we can 
conclude that the observers properly measure the ejected masses which are consistent 
with the super-Eddington theory describing the steady state after TNR. However, we 
theoreticians do not fully understand the mechanisms regulating the TNR threshold. 



which is obviously suppressed enough to allow more accumulation of mass before 
triggering a TNR. 

TABLE 1: The emerging "super-Eddington" picture. 

Phase Schematic Diagram (not to scale!) 

TNR. The nova eruption begins with a TNR. Initially, the 
system is dynamic. Some mass loss can occur from shocks or 
radiation driving. However, the amount of mass lost this way 
is insignificant compared to the total mass ejected. 

Formation of the super-Eddington steady state. Even before 
the wind reaches an equilibrium, a steady state is reached with 
shell burning, convection and a "porous" envelope. Before the 
wind builds up, any tenuous material already present above the 
critical point will accelerate to high velocities. 

Acceleration of a thick wind. Once the critical point is estab- 
lished in the envelope, a super-Eddington thick wind is accel- 
erated. With time, the steady state extends to include the wind. 
Since the wind is heavy, the radiation field is "used up" such 
that Ihe wind reaches slower ("primary") terminal velocilies. 

Stagnating wind conditions (Optional!). If the sonic radius 
shrinks fast enough, but the luminosity does not, conditions 
may form for a stagnating wind. That is, a wind too heavy to be 

pushed to r ^ oo is predicted such that there is no steady state. 
This could explain of the "oscillatory" type transition phase. 

Sub-Eddington state. Once a sub-Eddington state is re- 
established, the thick wind stops and the bare object emerges. 
The ejected gas is predominantly expanding with the "primary" 
velocity. A tenuous part has velocities typically twice as high, 
with a high dispersion - its outer parts accelerate to higher ve- 
locities. Thus, a Hubble flow wiU exist in the tenuous material. 
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